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The growth of large scale structure is studied in a universe containing both cold dark matter 
(CDM) and generalized Chaplygin gas (GCg). GCg is assumed to contribute only to the background 
evolution of the universe while the CDM component collapses and forms structures. We present 
some new analytical as well as numerical results for linear and non-linear growth in such model. The 
model passes the standard cosmological distance test without the need of a cosmological constant 
(LCDM) . But we find that the scenario is severely constrained by current observations of large scale 
structure. Any small deviations of the GCg parameters away from the standard Lambda dominated 
cosmology (LCDM) produces substantial suppression for the growth of structures. 



I. INTRODUCTION 



The current cosmological concordance model that is well supported by the latest microwave background and super- 
nova type la observations, includes both a large repulsive cosmological term and a cold dark matter (CDM) component. 
The nature and origin of these components is yet unknown, leaving room for novel theoretical explanations and new 
cosmological scenarios. In the standard cosmological scenario the cosmological term, A, is a true constant but, as is 
already suggested by the plethora of papers devoted to the subject, this need not to be the case. In addition, the 
difficulty of explaining the existence and value of the cosmological constant from fundamental physics, leaves plenty 
^ \ of room and motivation for phenomenological proposals. 

Instead of having a universe with a non-zero A, one can study cosmological models that, within the observational 
| boundaries, imitate the evolution of a universe with a cosmological constant and cold dark matter, the ACDM model. 
. Such models can have be the result of a modified Friedmann equation or contain a new type of matter with positive 
' energy density and negative pressure. Among these models are the large class of quintessence models (see e.g. 0), 
C*~) . the brane inspired models 0, Cardassian models and recently the Chaplygin^lJ gas models 0. 

The General Chaplygin gas (GCg) models are characterized by a perfect fluid with a non-standard equation of state 

o ■ p % 

^ \ where A is a positive (dimensionfull) constant. It is useful to define a dimensionless constant, A (also often denoted 
as A s ), by A = Pqq +ci 'A, where pc,o is the present day density of the CGg. 
. The original Chaplygin gas model Q corresponds to the choice a = 1. Such an equation of state can be the 
result of a scalar field with a non-standard kinetic term, e.g. the string theory motivated tachyon field 0, @- The 
generalization of the original Chaplygin gas, i.e. to cases where a ^ 1, was done in 0. Since the original Chaplygin 
gas scenario is reached as a special case of the general Chaplygin gas, we will concentrate on the latter. 

The cosmological behavior of the GCg is between a dark matter and a cosmological constant. At large densities, it 
behaves as pressureless dust where as at small densities, i.e. presently, it acts as a cosmological constant. As such, 
one can hope to have a unified dark matter candidate, i.e. one that can play both the role of cold dark matter and 
of the cosmological constant 0, . Taking a more conservative view, one can also consider a universe with both 
significant CDM and GCg components. The different cosmological models that include a (generalized) Chaplygin gas 
component can then be divided into two classes: models with and without a significant CDM component. Models 
that do not include a specific CDM component are often called unified dark matter (CDM) models. 

The unified dark matter models with a (generalized) Chaplygin gas have been studied in various works in view 
of cosmological observations. Cosmological perturbations of the GCg- fluid have been considered in 0,0]- 0- In 
|l2j it was claimed that the matter power spectrum strongly constrains the parameter space of unified dark matter 
models, essentially ruling them out. This result does not include the effect of the baryons, however, and is hence 
questionable 0| . Also, as it was pointed out in [l5l Il7|. the restriction due to the matter power spectrum may be 
avoided by e.g. allowing for entropy perturbations. To complicate the matter further, in [Tl ] it was shown that 
one must be careful in using linear theory to study the perturbations as linear theory may break down sooner than 
expected. 

Other observational measures of universes with a G C g- co mp onent include the SNIa constraints [1 [T3. fill IH. I20L 
IHIH, constraints coming from the CMB [1 HI EH 121121 E% and other astrophysical data [n llMl30f . 
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A significant fraction of the work so far has concentrated on the unified dark matter proposal. In this letter we 
concentrate on the more general case where there is in addition to the GCg-component, a significant cold dark matter 
fraction. The GCg plays the role of dark energy, from this point of view this model can be seen as a q uintessence 
model with a dynamical equation of state. GCg as dark energy has been studied in recent works El El El IM H3 • 
Interestingly, in |18l Il9| it was found by using SNIa data that a universe dominated with a Chaplygin gas was favored 
over a mixed model. On the other hand, in y3] it was concluded that a ACDM model is preferred. Also in j2|| , using 
WMAP data it was showed that a CDM model is disfavored and a mixed model, with a ACDM -like behavior, was 
preferred. The current bounds on the parameters are varied, depending on the type observational data. Using globular 
clusters to determine the age of the universe, in [3(| it was found that Qm > 0.2, A > 0.96 (a — 1 ), w here as for the 
lensing statistics the bounds are somewhat looser Qm < 0.45, ^4>0.72. Using the WMAP data, in [2|j the parameter 
bounds were found to be A > 0.8, < a < 0.2 (95%CL) and the original Chaplygin gas is ruled out. In [21|, X-rays 
from galaxy clusters, HST data and SNIa data were used to obtain A > 0.84, 0.273 < Q, M < 0.329 (95%C£). Using 
only SNIa data, it was found in El for the original Cg-model that A = 0.93t°;2o, < fl M < 0.35. 

In this work we study the mixed scenario, with CDM and GCg, from the point of view of large structure growth. 
We take a quintessence-like approach to the Chaplygin gas, i.e. only the CDM component is assumed to collapse 
and form structures. The Chaplygin gas acts only as dark energy and is visible, from the point of structure growth, 
only through the Friedmann equation. This is to be contrasted with the unified dark matter scenario where the Cg is 
also responsible for the growth of structures. The scenario considered here may be seen e.g. as a universe dominated 
by a scalar field with a non-standard kinetic term arising from string theory. Similar considerations for quintessence 
scenarios have been presented in [3l]], where large scale structure is used to test different quintessence models. 

By studying large scale structure growth from the point of view of gravitational collapse we probe the significance 
of the GCg dark energy component throughout the epoch during which large scale structures grow. The linear and 
non-linear growth can be related to observations and hence lead to observational constraints like galaxy number counts 
(see [32^. Our approach is similar to [HHil, where other non-standard cosmologies have been considered (also see 
|35| for related discussions). In Section 2 we review and study properties of GCg cosmology. The spherical collapse 
formalism is reviewed in Section 3 and applied to the GCg case, producing some analytical results. Numerical work 
is performed in Section 4. Discussion and conclusion is presented in Section 5. 



II. GENERALIZED CHAPLYGIN GAS COSMOLOGY 

Assuming in general that we have both a cold dark matter component and a GCg component, the Friedmann 
equation is 

H 2 ^Q 2 = K 2 (p M + PG ), (2) 

where k 2 — 8irG/3, pu is the energy density of cold dark matter and pa of the generalized Chaplygin gas. Note that 
there is no explicit cosmological constant and that we choose work in a flat universe. 
The continuity equation of the different components are 

p M + 3H(p M + Pm) = 
p G + 3H(p G +p G ) = 0, 

i.e. the two fluids interact only gravitationally. Using Eq. and the fact that in a matter dominated universe 
Pm = 0, we get two equations that can be solved straightforwardly: 

Pm,o 
Pm = — 3- 

PG = ^o(A+(l-A)a- 3 ( 1+Q )) T ^, (3) 

where the index refers to present day values (ao = 1) and A{= p G ^ +a ^ A) is a dimensionless constant, < A < 1. 
From Eq. J3J), we see that at early times, the GCg behaves as matter, p G ~ a -3 , whereas at late times it imitates a 
cosmological constant, pa ~ const. 

Using the expressions (01 and the Friedmann equation, Eq. @, one can look for solutions of a(t) or p(t). In the 
general case no analytical solution is found. As a special case, one sees that in a universe with only Chaplygin gas, 
i.e. pm = 0, the Friedmann equation is integrable and the solution is (up to a constant): 

2 1 + T ^ a 3( 1+Q ) yfa 1 1_ _± ^_ a 3( 1+Q)) (4) 



3\A+(l- A)a 3 ( 1 + Q )/ v 2(l + a)' 2(1 + a) ' 2(1 + a)' A— I 
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In the same special case, one can also solve for the energy density (again, up to a constant): 



l + 2a 



4a 



,l+£* 



1 3 Ak(1 + 2a) pl+ a 2Fl ( 2(1 + a) ' 2(1 + a) ' 



Using the solutions @, we rewrite the Fricdmann equation in terms of redshift as 

H 2 = Hl(n M {\ + zf + Q G (A+ (1 - A)(l + z) 3 ( 1+q ))t^) , 



(5) 



(6) 



where we have defined = k 2 Pm,o/Hq and ^g = k 2 pg : o/Hq. For a = we see that Eq. © reduces to the 
ACDM -model with f2 m = £Im + ^g(1 — A) and SI a = Vic A. Note that in a flat universe, our definitions imply that 
SijVf + f^G = 1- 

As a first test, we should compare distances in the GCg cosmology with the standard ACDM case. The luminosity 
distance is obtained from the (line-of-sight) comoving coordinate distance: r(z) = J dz' /H(z'). Fig. ^ shows 
predictions for different values of the model parameters for the corresponding apparent magnitude as compared to 
the SNIa results [28l |29| . It is apparent that current observations do not discriminate much between the different 
cosmological models. In other words: the GCg cosmologies seem to pass the standard cosmological distance test. It is 
then interesting to explore whether they also pass the observational constraints on the growth of linear and non-linear 
structures. 
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FIG. 1: Apparent magnitudes for SNIa for different cosmologies (Qm = 0.3, CIg = 0.7). 



III. GRAVITATIONAL COLLAPSE 



The gravitational evolution of an over-dense region is dependent on the background evolution. A rapidly expanding 
background will obviously slow down the collapse of region compared to a static background. The details of the 
collapse are obviously complicated in a general case so here we choose work in a shear-free approximation and 
consider spherically symmetric configurations, i.e. we use the spherical collapse model to study the growth of large 
scale structure. 

To study the dynamics of gravitational collapse, it is useful to define local density contrast as 



1. 



(7) 
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where p is the background energy density. Within the spherical colla pse approximation one can show that in a matter 
dominated universe (i.e. the collapsing matter has no pressure) [32, 33] 

where r\ — ln(a), H = H(p) is the background Hubble rate and H = H(S, p) is the local Hubble rate. 
Working within the framework of small perturbations we expand 8 as 



i=i i=i 



where 5q is the small perturbation. Similarly, we can also expand the H + H 2 -term in terms of 6 and then the whole 
RHS of Eq. © as 

(H + H 2 ) -{t+ H 2 )) = 3(1 + 6) J2 CnS n . (10) 



H 2 



n=l 

Using this and expanding the perturbation according to Eq. we find the linear equation, 

D'l + (2+^)D' 1 +3c 1 D 1 =0, (11) 
H 

the second order equation, 

D'i + (2 + ^)D' 2 - l(D[) 2 + 3 Cl D 2 + 6( Cl + c 2 )D\ = 0, (12) 
H 3 

and one can go on to arbitrary order by using the solutions to the lower order equations. 

The second order equation can be related to the skewness of the density field at large scales (see [34| for more 
details). For Gaussian initial conditions, 

S3 = 3^|. (13) 

In an Einstein-deSitter universe, 

cr = ~\ (14) 
d = 0, i = 2,3, ... (15) 

2 + l - \ <*> 

and the well known solution to the linear equation (|llfl is 

D x = Bie' 1 + B 2 e-i" 1 = B x a + B 2 a,-i. (17) 
The skewness can also be calculated analytically (see |34j): 

Si dS = ^ » 4.86. (18) 



A. Gravitational collapse in GCg 

The H + i? 2 -term in the GCg model is 



H + H 2 = - 



^ 2 (pa/ + Pg(1-3^L)) (19) 
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The gravitational dynamics in the GCg model obviously depend on the type of scenario i.e. whether we have a 
separate dark matter component or not and which components are fluctuating to form structures. 

Assuming that GCg only affects large scale structure growth through its effect on the background evolution, i.e. 
only the fluctuations of the CDM component are responsible for large scale structure, we can write pu = Pm(1 + S), 
Pg = Pg- Expanding according to (|10|l we find that 

ci=-\— — = 4 — (20) 

2 p M +PG 2 Q M + n G [\ + A(a 3 ( 1 + Q ) - 1)] — 

c 2 = (21) 

2 | = _n_|_3 ) = — (i-i — — i - - ) (22) 

H 2 2 K p2;{pm+Pg) 2 V Sl M a- 3 + Sl G [A + (1 - A)a~ 3 ( 1 + a )]^ J 

From Eqs (|20|l the small a limit, i.e. the early time limit, is easily read: 

ox = -\ ^ r- (23) 

2n M + n G (i-A)— 

2 + ^r = -. (24) 
H 2 2 v ' 

Comparing this to the linear equation in a ACDM -universe, 

„ l ^A f + 4a 3 A 3 Q M 

1 2 M + a 3 ft A 1 2 n M + a 3 il A ' 1 ' 

we see that a GCg- universe behaves fundamentally differently from the ACDM -universe: at early times, a <C 1, 
ACDM -universe reduces to the EdS-case, with solution given by Eq. ijTTI) . where as the GCg- universe does not. 
Accordingly, the solution to the linear equation in the GCg-case (in the small a limit) is 



D 1 =B ie i<-V^-^ + B2e i(-V^i-i)n i (2 6 ) 

where 

^ 

n M + n G (i-A)— 

Hence, EdS-case is recovered in the limits £l G = and A = 1. Compared to the EdS-case (and ACDM -universe), fluc- 
tuations therefore start to grow more slowly in a universe where GCg plays a role. Defining e = (1- Afl^ +a ^ G /n M , 
we see that the growing mode goes as (when e <C 1) 

£>i~a 1- * e . (28) 

On the other hand, in the large a limit, 

D» + ^M + ^oAf^jy _ 3 n M ^ = Qj 

2 f2 M + £l G A~a 3 2 Q M + fiG4^a 3 

which is the standard ACDM equation, Eq. 1)25(1 . with Q A = ^g^ 1+q • 
The second order equation reads at early times as 

D'i + \d' 2 - \iD 2 ^(D[f - 3£Dl = 0. (30) 

Substituting the growing mode from Eq. (|26|l , the solution to this equation is 



= 4 y i + 24$ - 1 - 21$ i (v ^-!^ 4(N /iT24e-i)^ + B jcyT+a^-i)^ (31) 

3 VI + 24$ - 1 - 18$ 



Hence, at early times, the skewness is given by 



S, = 3^ = 4^1-±^, (32) 
D? VTT24T- 1 - 18$ ' v ; 

which when fig = (i.e. $ = 1) reduces to the EdS-value, S^ ds = 34/7. At the other extremum point, SIm = or 
$ = 0, we reach S3 = 6. All the other values lie between these two limits. Numerical work (see later section) shows 
that Eq. (|32|l is a very good approximation until very late times, but even there the error is very small. 




B. Numerical Results 



In order to study the linear and non- linear evolution of the perturbation from early to late times, we must resort 
to numerical means. The linear and second order equations are solved numerically over the region of cosmological 
interest, from a = 10 -3 to a — 1 = a . The initial conditions are chosen such that in both cases, the growing part of the 
solution, Eqs Ij26|l and J2IJ, is followed at early times. This ensures that the behavior of the solutions is qualitatively 
similar to the ACDM -universe and that the ACDM -behavior is reached as A tends to unity. Normalization is chosen 
so that the amplitude of fluctuations is the same at present time, i.e. at a = 1. 

The linear factor is are plotted in Fig Note that the linear growth is depicted relative to a, i.e. to the growth 
rate in the EdS-universe (which is also equal to the growth rate in ACDM -universe at early times). In plotting the 
figures, we choose VLm = 0.3, which in a flat universe sets flc = 0.7. 

The linear growth rate confirms what was expected from the analytical results: linear growth occurs much more 
slowly than in a EdS (ACDM ) universe until at late times when the GCg acts as a cosmological constant that 
dominates, leading to even stronger suppression of growth. Note that the case A — 1.0 corresponds to the ACDM 
universe with S1m = 0.3, Ha = 0.7. The effect of changing a is much less significant as is shown in Fig. [3 (right panel). 

It is also interesting to study the time derivative of the above ratio: dF/dr = d(D/a)/dr 7 where dr is the comoving 
distance. This quantity provides a direct observable (see eg [3(| and references therein) through the integrated Sachs- 
Wolfe (ISW) effect (33,133. Numerical predictions for dF/dr are shown in Fig. [31 As expected deviations from the 
standard ACDM model are substantial for most of the relevant GCg parameters. 

The effect on the skewness is not as dramatic and is very well approximated by the analytical formula Eq. i)M2l) 

Further details on how these quantities can be compared to observations, can be found in |32l l33j and references 
therein. 



IV. DISCUSSION 



We have considered a simple variation on the Generalized Chaplygin gas cosmology, where in addition to the GCg- 
component there is a significant cold dark matter fraction. The GCg plays the role of dark energy, but not in a unified 
way (see Introduction) meaning that the GCg does not collapse gravitationally. At first look this model seem viable, 
as it has quite a standard looking Friedmann equation, see Eq. (|2J) , and cosmological distances, see Fig. ^ However, 
this is not the case for the growth of structures. 

As we have seen both analytically and numerically, the growth of large scale structure is fundamentally different 
from the ACDM scenario due to the properties of the Generalized Chaplygin gas. Unlike in the ACDM -model, 
where at early times the equation governing linear growth equals that of the EdS-case, the GCg behaves differently 
at early times. From Eq. we see that for A ^ 1, pa : o 0, linear growth will be suppressed relative to the EdS 
scenario, and to a ACDM universe which follows the EdS-solution until late times. We can see in Fig. [3 that any 
small deviations of the GCg parameters from the ACDM values (a = 1 and ^4=1) produced substantial suppression 
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FIG. 3: Derivative of the linear growth: F = d(D/a)/dr for a = 1.0, A = 1,0.95,0.9,0.7,0.5 (left panel), A = 0.9, a = 
0,0.2, 1.0(right panel) 

of linear growth. For example, even the case A = 0.99 produces a much higher amplitude of fluctuations at decoupling 
(a ~ 10 -3 ) than the ACDM cosmology (when normalized at a = 1). Such differences in normalization are hard to 
reconcile with our current understanding of CMB anisotropies (see eg |32l l40j). The final numbers dependent on the 
adopted model for CMB transfer function and matter content, but it is already apparent that we will need to tune 
the GCg parameters very close to the ACDM model. As an independent test we have also presented the predictions 
of the GCg model for the ISW effect, see FigEland the skewness of the density field, see Eq.(jj2J). Again, the ISW 
effect is much larger even for A — 0.99 compared to a ACDM universe. It should be stressed that in the scenario 
considered here, the GCg only acts through the background evolution and does not undergo gravitational collapse. 
Hence, these parameter constraints do not apply to the unified dark matter models with a GCg component. 

Observations seem to favor values of the GCg parameter which makes it equivalent in practice to a cosmological 
constant. We conclude therefore that using the GCg as a form of dark energy does not seem to provide any apparent 
advantage over the more simple ACDM model, which has fewer parameters. It is nevertheless clear that future 
observations will be able to provide tighter constraints on the quantities we have studied here. If the ACDM model 
fails to match such observations, the GCg models could be an useful alternative to consider. 
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